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Abstract. We provide a brief review of thermohaline physics and why it is a candidate ex- 
tra mixing mechanism during the red giant branch (RGB). We discuss how thermohaline 
mixing (also called S/-i mixing) during the RGB due to 3 He burning, is more complicated 
than the operation of thermohaline mixing in other stellar contexts (such as following accre- 
tion from a binary companion). We try to use observations of carbon depletion in globular 
clusters to help constrain the formalism and the diffusion coefficient or mixing velocity that 
should be used in stellar models. We are able to match the spread of carbon depletion for 
metal poor field giants but are unable to do so for cluster giants, which may show evidence 
of mixing prior to even the first dredge-up event. 
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1. Introduction 

The need for extra mixing on the RGB is ob- 
servationally well established. Any mechanism 
(or the combined effect of multiple mecha- 
nisms) must meet the following requirements: 



L. 



It must occur after the luminosity bump 
and continue to operat e until near the 
tip of the RGB. dGilroy & Brown} 

120041 
iMartell et all 



1991; 

Smith jfcMartel 



2008) 



Weiss & Charbonnel, 
2003T~ 



2. It must occur over a range of mass es and 
metallicities. (Smilian ic et all l2009h 

3. It must deplete 7 Li. dCharbonnel et all 
1 1 9981: LSmilianic et al l 120091) ' 

4. Itmust deplete 3 He. dDearborn et"aiT,ll986t 
lHata et aUll995tlDearborn et all 1 19961) 



It must lower the 12 C/ 13 C ratio 
(ICharbonnelill994ll996l) 
It must deplete the carbon abundance 
and increase th e nit r ogen abundance. 
( Smi lianic et aD. 120091: 



2003; Martell et ah, 2008) 



Smith & Martell, 
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These criteria suggest that, in order for 
theory to remain consistent with observations, 
material must be mixed through radiative re- 
gions, processed by the H-shell, and mixed 
back into the enevelope. This requirement is 
often referred to as deep mixing because mix- 
ing deeper than the formal convective bound- 
ary into the radiative zones will lead to ma- 
terial being exposed to regions of higher tem- 
perature and will result in the required addi- 
tional processing. In general the 12 C/ 13 C ratio 
is used t o probe the effici e ncy of first dredge up 
(FDU); dDearborn et all Il975b iTomkin etall 
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119761) and is also used as a tracer of the extent 
of deep mixing. A good example of this was 
Sweigart & MengeH d 19791) who were the first 
to use the isotopes to investigate the role of 
rotational mix i ng on the RGB. More recently 
iPalacios et al.l {2006) have shown that, whilst 
rotation does reduce the 12 C/ 13 C isotopic ra- 
tio, it is unable to explain the values seen in 
giant photospheres. Although it is understood 
that extra mixing must take place, only recently 
has a mechanism (i.e. thermohaline mixing) 
been discovered that can poten tially satisfy all 
of the aforementioned criteria (Eg gleton et al.L 
120061) . 

2. Thermohaline Mixing In Stars 

2.1. Historical Overview 

Thermohaline mixin g was firs t studied in the 
Earth's oceans bv LSternl dl960h where strati- 
fied warm salty water sits upon a cool unsalted 
layer. The layers are initially stable. However, 
heat diffuses more quickly than composition 
so the warmer layers cool. Now they are sim- 
ply denser than the material underneath, and 
a turnover is initiated via the formation of 
lengthy "fingers" of cooler salty water reach- 
ing down into the cold fresh water. This dis- 
places cool fresh water upwards, and a mixing 
occurs. On a slower timescale the salt diffuses 
out of the salty cool water to reach a new salti- 
ness in the mixed region. 

This double diffusive m ixing was first ap- 
plied t o a stell a r cont ext by IStothers & Simonl 
dl969b . |Ulrichl dl972l) applied th is to a perfect 
gas and Kiorjenhahn et al extended this 

to allow for a non-perfect gas which included 
radiation pressure and degeneracy. There were 
two obvious situations in which they applied 
thermohaline mixing. Firstly, during pre-main 
sequence contraction when in-situ 3 He burning 
lowers the local mean molecular weight, // be- 
cause the reaction 



3 He ( 3 He, 2p) 4 He 



(1) 



molecular weight. This was found to have lit- 
tle effect, due to the short pre-main-sequence 
time scale and the fact the star becomes fully 
convective before reaching the Zero Age Main 
Sequence (ZAMS). The second case consis- 
dered was during the core flash, when during 
off-centre He ignition, carbon-rich material sits 
upon helium-rich material. This also was con- 
sidered to have little effect on the evolution 
primarly due to the uncertainty of competing 
timescales. The mixing must occur before the 
star settles down to q uiescent helium burning. 
Eggle ton et al.l (2006) also showed that a small 
amount of overhooting inwards could remove 
the molecular weight inversion on a dynamical 
timescale. 



2.2. Application to the RGB 



produces more particles than it destroys. The 
mixing is determined by the competition of 
the heat diffusion and the difference in com- 
position but it is driven by the change in local 



lEggleton etail d2006l) used a 3D hydrody- 
namical stellar code dDearborn et all 120061) 
to show an instability develops due to 3 He 
burning along the RGB, an instability that 
ICharbonnel & Zahnl ( 20071) model as thermo- 
haline mixing. Charbon nel & Zahnl (2007) and 
Eggle ton etafl d2008l) have modelled this 3 He 
burning driven instability in their ID codes and 
demonstrated the significant effect the mix- 
ing has during the RGB. Following FDU, the 
convective envelope recedes, leaving behind 
a homogeneous region. Any composition and 
molecular weight gradient has been removed 
due to the convective mixing. As the hydrogen 
burning shell begins to advance, 3 He begins to 
burn. From Equation Q] it can be seen that this 
reacti on creates a loc a l mo lecular weight inver- 
sion; Eggle ton et ail d2008l) found its magni- 
tude to be of the order Afi/fi ~ 10~ 4 . Although 
the inversion seems small, convection is in fact 
driven by a similarly small superadiabaticity. 
Usually such a small change in the local molec- 
ular weight would have almost no effect, as it 
would be swamped by the existing yu gradient 
produced by the burning of other species. It is 
this unique situation where 3 He begins to burn 
before the other species and the fact that first 
dredge up has homogenised the region that al- 
lows the inversion to develop. 

Although there is no salt in the star the pro- 
cess is doubly diffusive and thus labelled ther- 
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mohaline mixing. The authors refer to this as 
6fi mixing to emphasise that the mechanism 
that drives the mixing and the fact it is more 
complex than the other examples of thermo- 
haline mixing. As 3 He burns, a parcel forms 
that is hotter and has lower molecular weight 
than its surroundings. It quickly expands (and 
begins to cool) in order to establish pressure 
equilibrium. The expansion reduces the density 
and therefore the element becomes buoyant. 
The parcel rises until it finds an equilibrium 
point where the external pressure and density 
are equal to that inside the bubble. This is ex- 
pected to be a small displacement which occurs 
on a dynamical timescale. 

As the molecular weight inside the bubble 
is lower than its surroundings the equilibrium 
point must correspond to a place where the ex- 
ternal temperature is higher than that of the 
bubble. The temperature inside the bubble will 
be lower than its surroundings: 



Ti 
T 



Ho' 



(2) 



where subscript i denotes the inside of the bub- 
ble and subscript o denotes the surroundings. 
As heat begins to diffuse into the parcel, we 
expect layers will start to strip off in the form 
of long fingers. It is this secondary mixing 
that governs the overall mixing timescale. The 
mixing cycles in fresh 3 He from the envelope 
reservoir, while CN-processed material is 
cyc led into the con ve ction zone. 

lEggleton et all d2008i EDL hereafter) 
found that this mixing satisfies the criteria 
outlined in Section 1. The level of depletion 
of the carbon isotopes will depend on the 
efficiency of the mechanism. EDL estimated 
the mixing speed and with their formula 
for the diffusion coefficient found that a 
window of three orders of magnitude in the 
mixing velocity can lead to o bserved levels 
of 12 C/ 13 C and 3 He depletion. lUlrichl (Il972l) 
and iKippenhahn. Ruschenplatt. & Thomas! 
(1980) use essentially the same formula for 
the diffusion coefficient (UKRT hereafter) but 
their geometric c oefficients vary by tw o orders 
of magnitude. Charbonnel & Z ahnl d2007l) 
have applied the UKRT mixi ng to the RGB. 
Both lCharbonnel & Zahnl d2007h and EDL see 



the 12 C/ 13 C ratio is reduced to similar levels. 
In this study we will attempt to use globular 
cluster observations to constrain both the form 
of the diffusion coefficient and the mixing 
velocity. 



3. The Mixing Speed 

In order to implement 5/u mixing into our ID 
codes we must consider the following: 

1 . Which formalism should be used? Here we 
will limit our investigation to the EDL and 
UKRT prescriptions for the diffusion coef- 
ficient. 

2. Once the preferred formalism is identified 
what mixing velocity is needed to match 
observations? What values do we use for 
any free parameters? 

3. The 12 C/ 13 C ratio is generally used as a 
proxy to probe the extent of mixing. This 
quickly saturates in low metallicity stars 
and therefore could be misleading. Is there 
a better way to try to constrain the velocity? 

EDL postulated the following formula 
based on the velocities from their 3D code in 
analogy with the existing convective formalism 
in their code: 



( (fi - fXain) if (k > k mia ) 
D = Wk; " (3) 







if (k < k min ), 



where yu m i n is the smallest value of ji in 
the current model, k the mesh point number, 
counted outwards from the centre, r is the ra- 
dial coordinate, F; nv is a constant which is se- 
lected to obtain the desired mixing efficiency 
and fnuciear is an estimate of the nuclear evolu- 
tion timescale (see EDL). 

This formulation ensured the correct re- 
gion was mixed but ensures the mixing is for- 
mally zero at the position where fi has its min- 
imum even though it should presumably be the 
most efficient at this point. EDL give upper 
and lower estimates for the mixing velocity and 
find that they can alter the speed by three orders 
of magnitude and still produce the observed 
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levels of 12 C/ 13 C and 3 He. lCharbonnel & Zahnl 
d2007l) adopt the UKRT formula 

D '- C,K ^)^) f ° r ^ <0 ' W 

where V = (dlnT/dlnP), <p 
(d\np/d\nfi) P j, 6 = -(dlnp/dlnT)^ 
and C, is the geometric factor. 

Empirical studies of flu ids in laboratory 
conditions led Ulrichj (1 19721) to determine that 
C t ~ 1000. He saw the development of long 
salt fingers with lengths that were larger than 
their diameters, which led to efficient mixing. 
Kippenhahn on the other hand envisaged the 
classical picture where mixing is due to blobs 
and thus determined C t ~ 10. 

We have run stellar models of various 
masses, with both EDL and UKRT mixing. We 
tested different values of F lm and C t in order to 
alter the efficiency of mixing. To test our mod- 
els for the extra mixing we chose to use the 
carbon abu ndance as a func t ion of My as deter- 
mined bv ISmith & Martelll d2003l) . They plot- 
ted carbon abundance as a function of visual 
magnitude for a variety of globular clusters. 
In doing so they were able to clearly demon- 
strate the depletion of carbon along the RGB. 
Globular cluster have always been an excel- 
lent test bed for stellar theory and by trying to 
match the carbon depletion for various red gi- 
ant branches we have an alternative abundance 
test for mixing efficiency. 

4. Results 

In Figure 1 we plot carbon abundance for stars 
in the Galactic glob ular cluster M92 and th e 
Galactic halo from ISmith & Martelll d2003l) . 
Open circles denote galactic field giants whose 
metallicity rang es from -2.0 < fFe/HI < -1.0 
also taken from ISmith & Martelll (120031) . The 
filled circles correspond to RGB stars in M92. 
In both the field and the halo it is immedi- 
ately obvious that there is carbon depletion as 
stars ascend the giant branch. If our models are 
able to match the carbon depletion we may be 
able to constrain the thermohaline mixing for- 
malism and velocity. Another thing to notice 
before turning to the models is the spread in 



carbon for a given visual magnitude. We at- 
tribute this to the primordial abundances of the 
cluster with the most C-rich at a given mag- 
nitude being the "normal" stars. The spread 
in C at a given magnitude is assumed to be 
of primordial origin as is the case with many 
other globular clusters. Our primary aim is 
to match the level of carbon depletion. That 
is, we are concerned with matching the de- 
crease in the upper and lower limits of the 
[C/Fe] values, as a function of magnitude. The 
solid and da shed lines were computed using 
MONSTAR (ICampbell & Lattanziol200"8h . We 
have evolved a 0.8 M Q and a 0.9 M star until 
the core flash. These masses straddle the age 
limits of stars in this cluster. A metallicity of 
Z=0.0001 was used to match that of the M92 
where [Fe/Hl^-2.2 dBellman et alll200ll) . The 
EDL mixing quickly destroys the 3 He without 
significantly altering the FDU values of car- 
bon. We believe this model is not mixing to 
high enough temperatures. As mentioned in the 
previous section, the mixing speed is formally 
zero at the position where fi has its minimum. 
By not mixing at the minimum properly the yu 
profile is being affected and carbon is not be- 
ing exposed to the required temperature in ei- 
ther model. 

The dot-dashed line is a model computed 
using the Eggleton code, dEggletonl 11971b 
IStancliffe & Eldridgei 120091) . It too is of mass 
0.8 M and corresponds to the metallicity of 
M92 however it is run without mass loss. 
Running without mass loss here will result in 
less carbon depletion than we would otherwise 
expect, therefore will serve as a lower limit for 
the depletion of carbon. An EDL style formula 
for the diffusion coefficient is used in this cal- 
culation, that is there is a dependence on the 
position where yu reaches it's minimum. The 
Fi nv here was calibrated such that a 1.5 M , 
Z=0.0001 model gave the same level of carbon 
depletion on the RGB as a 1.5 M Z=0.0001 
mo del with UKRT mixing where C t =1000, 
(see Stancliffe l2010l for more detail). 

The dotted line is a model with a UKRT 
prescription taken from lStancliffe et al.l (2009) 
where C t =1000. This also was run without 
mass loss. The UKRT mixing is a local formal- 
ism that is dependent on the pi gradient. Unlike 
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o o Halo Field Giants 

• • Red Giants in M92 

M08Z0001 Finv=100 

.... M09Z0001 Finv=100 

RJSm08z0001 (EDL like) 

RJSm08z0001 (Kipp) 



Mv 



Fig. 1. We plot the carbon abundance [C/Fe] as a function of M v for observed giants and our models. We 
include cluster giants from M92 (solid circles) along with Galactic halo field giants (ope n circles) where 
the me tallicity covers [Fe/H] =-1 to -2. Both sets of observational data were taken front Smith & Martell 
( 2003). The solid line and the bashed line correspond to models run with MONSTAR. The solid black line 
is the evolutionary track for a 0.8 M Q star up until the core flash. The dashed line is the evolution of a 0.9 
M o star. The dotted and dot-d ashed lines are models run with the Eggleton code which has been modified 
bv lstancliffe & Eldridgd O>009). The dot-dashed line corresponds to the evolution of a 0.8 M Q star using an 
EDL formalism. The dotted line is for a 0.8 M Q star using a UKRT style mixing and C t = 1000. 



in the EDL case this translates to the mixing 
being more efficient at the position where fi 
reaches its minimum. In both cases carbon is 
brought down from the envelope but here it 
is mixed to the position of lowest molecular 
weight and hence exposed to the shell much 
faster. The high temperature gradient ensures 
that mixing only a little deeper will see the car- 
bon undergo larger depletion. This is of course 
all dependent on the amount of 3 He available to 
drive the mixing. We see that the UKRT mix- 
ing can lead to levels of depletion seen in the 
field giants. Given that the field giants and M92 
stars are of similar age and metallicity it is in- 
teresting the cluster stars undergo more sub- 
stantial depletion. We defer the discussion of 
why this is to subsequent work. 



5. Conclusion 

Our initial motivation behind this paper was 
to use the observed variation of carbon abun- 
dances on the giant branch to help constrain 
some of the uncertainties present in the ther- 
mohaline mixing which we believe is operating 
during the red-giant phase. Drawn by the best 
data being available for M92, we chose this as 
our first attempt to fit the observations. The fact 
that we have failed in our aim has nevertheless 
taught us three important things: 

1. The functional form of the diffusion co- 
efficient strongly influences the depletion 
of carbon. 

2. Comparing the carbon isotope ratio is not 
necessarily useful because it saturates at 
the equilibrium value of about four while 
C continues to burn into N. 

3. The carbon abundances in M92 may pro- 
vide a very serious challenge for stellar 
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evolution, independent of any deep-mixing 
mechanism. At the same time it could in 
fact be telling us something very important 
about the deep mixing process. 

Concerning the first point we note that the sim- 
ple formula used by EDL causes an initially 
rapid depletion and then a leveling off, which 
does not seem to match the observations for 
metal-poor globular clusters. The UKRT de- 
scription results in a more gradual depletion 
and may be a better description. Neither de- 
pletes the carbon by enough to match the ob- 
servations, but we note that the Eggleton mod- 
els here are run without mass loss which will 
exaggerate the discrepancy. 

We believe that the third point is more fun- 
damental. The data for M92 clearly show de- 
pletion in [C/Fe] for stars with magnitudes My 
> 1. Note that standard stellar evolution pre- 
dicts that the first dredge-up does not produce 
observable abundance changes for these stars, 
and that this dredge-up does not finish until a 
magnitude My ~ +0.5. By this stage in the evo- 
lution, the stars are already showing depletions 
of C of order 0.5 dex. Further, the bump in the 
luminosity function (hereafter LF b ump) is ob- 
serve d to be at My ~ -0.4 dFusi Pecci et all 
1990). According to the usual ideas, deep- 
mixing (by whatever the mechanism) is inhib- 
ited until the star reaches the LF bump and 
the advancing H-shell removes the molecu- 
lar weight discontinuity left behind by the re- 
ceding convective envelope at the end of first 
dredge-up. In the case of M92 the stars on the 
giant branch have already depleted their [C/Fe] 
by about 0.8 dex when they reach this stage. If 
we have to postulate that some form of mix- 
ing begins sufficiently early to produce this de- 
pletion, then the mixing must necessarily re- 
move the abundance discontinuity that is itself 
responsible for the observed LF bump! The re- 
sulting contradiction produces, in our view, a 
serious problem for stellar astrophysics. 

It is worth noting that the LF bump in M92 
is not as cle arly visible as it i s in m ore metal- 
rich clusters. iFusi Pecci et al.l (11990 ) had to co- 
add data for three very similar clusters to make 
it visible in the d ata. Indeed, recent work by 
iPaust et all d2007l) provides little evidence for 



a bump in the observed LF of M92. These au- 
thors show that even the t heoretically pre dicted 
bump is small (see also ISweigartlll978l) . We 
are left trying to identify cause and effect: is 
the reduced bump the result of a reduced dis- 
continuity in the molecular weight in this case, 
which is not enough to prevent mixing before 
the disconintuity is erased by nuclear burning? 
Or does some mixing begin before the bump is 
reached, with the necessity of that mixing re- 
ducing the molecular weight discontinuity? 

We note that we are not the first to have no- 
ticed this problem, as it has be e n discussed by 
fat lea st) iMartell et alJ d2008h. iBellman et all 
(1200 ll) . and lLanger et al I dl986h . However, the 
data in Figure 1 are compiled from various 
sources and this presents a uniformity problem. 
Offsets by 0.3 dex are possible (G Smith, pri- 
vate communication) and could be the cause 
of the apparent contradiction. Certainly to use 
M92 as a constraint for dfi mixing requires 
a homogeneous set of data covering a wide 
range of luminosities. Such data are simply not 
available at present, but would prove extremely 
valuable. 
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